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Abstract 

We report on a study of direct CP violation in the decay D'^*^^K^ with a Dalitz 

analysis of the D'^ Kstt^tt'^ decay using a sample of 227 million BB pairs collected by 
the BABAR detector. Reference to the charge-conjugate state is implied here. We constrain 

the amplitude ratio rB= \'a(b-^d"k-]\ ^ ^-^^ ^^^'^ confidence level and r^ = 

[^(B-Z^i^tox- j{ = 0-155^0 077 ^ 0.040 ± 0.020 and we measure the relative strong phase Sb = 

(114 ± 41 ± 8 ± 10)° and 6*g = (303 ± 34 ± 14 ± 10)° between the amplitudes A{B- d''*^°K- ) 
and A{B- ^ D^*'>°K-). From these samples we measure 7 = (70 ± 26 ± 10 ± 10)°. The first 
error is statistical, the second error accounts for experimental uncertainties and the third error 
reflects the Dalitz model uncertainty. For this preliminary result we have quoted confidence 
intervals obtained with a Bayesian technique assuming a uniform prior in rs, 7 and 5b- 
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1 INTRODUCTION 



In the past years CP violation in the B meson system has been clearly established ^ and although 
there is good agreement with the expectations of the Standard Model, further measurements of 
CP violation in B decays are needed to over-constrain the unitarity triangle and look for New 
Physics effects. A crucial test will be represented by the measurement of 7, which is the complex 
phase of the Cabibbo-Kobayashi-Maskawa |2j quark mixing matrix element V^f, in the Wolfenstein 
parameterization [HI. 

Various methods using B D^K^ decays ||4j have been proposed to measure the unitarity 
triangle angle 7, all exploiting the fact that a B~ can decay into a D^K~ final state via a 6 — > c 
transition or into a final state via a b ^ u transition. CP violation can be detected if the 

and D decay into the same final state. The measurement of direct CP violation is sensitive 
to the phase difference between Vub and Vcb and thus to the angle 7. Most of the experimental 
methods to extract 7 can be grouped in two categories: the and D decay into a CP eigenstate 
[H]; or the decays to a common flavor state, either through a Cabibbo-allowed or a doubly 
Cabibbo-suppressed mode [Hj. The measurement of 7 in both methods also requires the knowledge 
of r^, the magnitude of the ratio of the amplitudes A{B~ D K~) and A{B~ D^K~) and of 
their relative strong phase 5b, which can be obtained from data. 

In this paper we report on a measurement of direct CP violation in B^ D^*^^K^ based on the 
analysis of the Dalitz distribution of the three-body decay i^svr^vr"'" j7j. The advantage of this 

method is that it involves the entire resonant substructure of the three-body decay, with Cabibbo- 
allowed and doubly Cabibbo-suppressed amplitudes interfering directly. It is therefore expected to 
have a higher statistical precision than the methods outlined above. Results of an analysis based 
on this procedure were reported by the Belle Collaboration in jH|. From the combination of the 
B~ D^K- and B' D*^K- mode they obtain the value 7 = 77°l|^o ±13°±11° where the flrst 
error is statistical, the second is experimental systematics and the third is model uncertainty. They 
also obtain a value of rs= 0.26lg:J^± 0.03 ±0.04 for ^ D'^R- and r*^ = 0.20l[]:J?± 0.02 ±0.04 
fovB~^D*^K-. 

1.1 Analysis outline 

The B and B^ decay amplitudes for the B^ D^*^'^K^ and Kst^^t^^ decays can be 

written assuming no CP asymmetry in D decays as : 

M_{m'i,m\) = \A{B- ^ D^K-)\\f(m'i,m\)+rBe'^^''-^'^ f{m\,m'i)], (1) 

M+{m^_,m\) = \A{B+^D^K+)\\f{m\,m^_)+rBe'^^^^^^f{m^_,m\)\, (2) 

where m?_ and are the squared invariant masses of the Kstt" and Kstt^ combinations respec- 
tively and /(m^,m^) is the amplitude of the KsTr~ir~^ decay. 

Given a known /, the bi-dimensional Dalitz (rn-?_,m^) distributions for B~ and B~^ can be 
simultaneously fitted to |M_(m?_, m^)p and |M+(rn-^, m^)p respectively. A maximum likelihood 
technique may be used to estimate r^, 6b, and 7. Since the measurement of 7 arises from the 
interference in Eq.^and Eq.[21 the uncertainty in the knowledge of the complex form of f{rn?L,rn\) 
can lead to a systematic uncertainty. A model describing the — > Kstt~t^~^ decay in terms of 
two-body amplitudes has been assumed in this analysis. This model has been characterized using 
a high statistics flavor tagged sample (Z)*+ — > D^irf ), obtained from e~^e~ — > cc events as 
described in Section IH 
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A similar analysis is also performed using B~ — > D*"K~ decays, and 7 is extracted along with 
the amplitude ratio r^j and strong phase difference 6^^ taking into account the effective strong 
phase shift of n radians between the D*^ D^n^ and D*^ D^j channels By convention 5^ 
is the strong phase of D*^ D^ir^ decay mode. 

2 THE BABAR DETECTOR AND DATASET 

The analysis is based on a sample of 227 million BB pairs collected by the BABAR detector at 
the SLAC PEP-II e^e~ asymmetric-energy storage ring. BABAR is a solenoidal detector optimized 
for the asymmetric-energy beams at PEP-II and is described in ^|. We summarize briefly the 
components that are crucial to this analysis. Charged-particle tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). In addition to providing precise 
spatial hits for tracking, the SVT and DCH also measure the ionization energy loss (dE/dx), which 
is used for particle identification of low-momentum charged particles. At higher momenta (p > 0.7 
GeV/c) pions and kaons are identified by Cherenkov radiation detected in a ring-imaging device 
(DIRC). The typical separation between pions and kaons varies from 8a at 2 GeV/c to 2.5cj at 
4 GeV/c. Neutral cluster (photon) positions and energies are measured with an electromagnetic 
calorimeter (EMC) consisting of 6580 thallium-doped Csl crystals. Candidate vr^ mesons are re- 
constructed as pairs of photons, spatially separated in the EMC, with an invariant mass within 3a 
of the vr'^ mass. These systems are mounted inside a 1.5-T solenoidal super-conducting magnet. 



3 EVENT SELECTION 

We reconstruct the decays B" D^K~ and B~ D*^K- with D*^ D^tt^ , D^j. A larger 
sample of B~ — > D^*^'^Tr~ is also reconstructed and is used as a control sample to determine the 
Probability Density Function (PDF) of the discriminating variables used in the likelihood fit for 7. 

candidates are reconstructed in the ET^Tr^vr" final state with the Ks reconstructed from pairs 
of oppositely charged pions with an invariant mass within 9 MeV/c^ of the nominal Ks mass 
The two pions are constrained to originate from the same point. The angle between the Ks 
line of flight and its momentum is required to satisfy the condition cos axg >0.99. D° candidates 
are selected by making all possible combinations of the Ks candidate and two oppositely charged 
pions with an invariant mass within 12 MeV/c^ of the nominal mass. 

The photon candidates for D*^ — > D^j are reconstructed from clusters in the electromagnetic 
calorimeter with energy greater than 30 MeV and consistent with a photon shower proflle. We 
select candidates from pairs of photon candidates and require 115 < m(77) < 150 MeV/c^ and 
with total energy greater than 70 MeV. To improve the momentum resolution, the vr*^ candidates 
are kinematically fltted with their mass constrained to the nominal vr*^ mass. The candidates 
are combined with a low energy vr^ or 7. The D*^-D^ mass difference Am is required to be within 
2.5 (10) MeV/c^ of the nominal Am for D^-k^{'^). 

A B~ candidate is obtained by combining a D^*^'^ candidate with a track ("bachelor" track) 
identified as a kaon as described in jlJlj. We improve the momentum resolution for the daughters 
by applying a kinematic mass constraint. For every B candidate two standard variables are defined. 
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the beam-energy-substituted mass rngs = y {\s + pq ■ ps)"^ / Eq — p\ and the energy difference 
A£^ = E'^ — ^\/s, where the asterisk denotes the CM frame, s is the square of the total energy 
in the CM frame, p and E are, respectively, momentum and energy, and the subscripts and B 
refer to T{AS) and B^, respectively. The resolutions, evaluated on simulated signal events, are 
2.6 MeV/c^ and 17 MeV for mES and /S.E, respectively. 
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Figure 1: B' D^R- (top left), B' D*^[D^tt^)K^ (top right) and B' D*°{D^-f)K- 
(bottom) niES distribution in the AE region [—30,30] MeV in the sample of 211 million BB pairs. 
The signal contribution is shown in red, B~ ^ n*0^- 
in magenta. 



D TT in blue, generic BB in green, and continuum 



To distinguish between BB and continuum events the following topological variables are used: 
cos^Ij, where 0^ is the polar angle of the B candidate with respect to the beam axis in the CM 
frame; Lq = J2iPi ^-nd L2 = J2iPi cos^ 9i calculated in the CM frame, where pi and 9i are the mo- 
menta and the angles of tracks and neutral clusters not used to reconstruct the B candidate with 



10 



respect to its thrust axis; cos 9t, where 9t is the angle in the CM frame between the thrust axes of 
the B candidate and of the rest of the event. B candidates are selected by requiring | cos 6t\< 0.8. 
Under this condition a Fisher discriminant is constructed from the variables discussed above. 
This Fisher discriminant is used in the likelihood fit to help distinguish between signal and con- 
tinuum e^e~ qq {q = u,d,s,c) events. The PDF of J-' is obtained for continuum events using 
niES sideband data. The B~ D^*^^Tr~ events are used to obtain the PDF for BB and signal 
events. 

Finally, the B~ D^*^'^K~ sample events are required to satisfy rriEs > 5.2 GeV/c^, \AE\ < 30 
MeV. The D^*'i^'n:- candidates are selected using criteria similar to those applied for 

B- D^*^^K~ but requiring that the bachelor pion not be consistent with the kaon hypothe- 
sis. The overall reconstruction efficiencies are 18%, 4.3%, 8.1% for the D^K~ , D*^(D^it^)K~ , and 
D*^{D^j)K~ decay modes, respectively. Fig. ^shows the rriES distributions after all the selection 
criteria are applied. 

4 DETERMINATION OF A'^tt 7r+ DECAY MODEL 

The amplitude /(m^, m^) has been constructed from a Dalitz analysis of a 97% pure flavor tagged 
sample obtained from 81496 D*^ D^n'^ events corresponding to a luminosity of 91.5 fb~^ 
(Fig. 12)). The Dalitz (m^,m^) distribution (Fig. E)) is fitted in the context of the isobar formalism 
described in In this formalism the amplitude / can be written as a sum of two-body decay 
matrix elements and a non-resonant term according to the expression 

/ = On^e*'^"'- + Urare"''^As{Ks7r~7:^\r). (3) 

Each term of the sum is parameterized with an amplitude and a phase. The factor As{KsTT~iT~^\r) 
gives the Lorentz invariant expression for the matrix element of a meson decaying into Kstt^tt^ 
through an intermediate resonance r as a function of the position in the Dalitz plot. It is, in general, 
parameterized by a relativistic Breit-Wigner with a functional form dependent on the spin of the 
resonance. For the p a more complex parametrization is used as suggested in [13]. We fit the Dalitz 
distribution with a model consisting of 13 resonances leading to 17 two-body decay amplitudes and 
phases (Table Of the 13 resonances eight involve a Ks plus a vrvr resonance and the remain- 
ing five are made of a {Kstt~) resonance plus a vr"*". We also include the corresponding doubly 
Cabibbo-suppressed amplitudes for most of the {Kstt~) vr"^ decays. All the resonances considered 
in this model are well established except for the two scalar tttt resonances, ai and o"2, whose masses 
and widths are obtained from our sample. Those are introduced in order to obtain an acceptable 
fit to the data, but their existence as true, scalar particles is a matter outside the scope of the paper. 

An unbinned maximum likelihood fit is performed to measure the amplitudes anr,CLr and the 
phases (pnr , (l>r ■ The fit fraction for each decay channel is defined as the integral of a single com- 
ponent divided by the coherent sum of all components. The results of the fit are shown in Fig. |31 
Amplitudes, phases and fit fractions as obtained by the likelihood fit are reported in Tabled 

We estimate the goodness of the fit for our model with a fit using adaptive binning of the 
Dalitz plot. We obtain xVd-o.f. = 3824/(3054-32) = 1.27. 

To illustrate the region of the Dalitz plot most sensitive to 7 measurement, we show in Fig. |1] 
the distribution of simulated B~ — > D^K~ events based on our Dalitz model, where each event is 
given a weight of ^-J where C is the likelihood function described in the following section. The 
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Resonance 


Amplitude 


Phase 


Fraction 


Mass 


Width 


Functional 






(degrees) 


(%) 


MeV/c2 


MeV/c2 


form 


K*{892) 


1.777 ±0.018 


131.0 ±0.81 


58.51 


891.66 


50.8 


BW 


p'J(770) 


1 (fixed) 


(fixed) 


22.33 


775.8 


146.4 


GS 


i^*(892) DCS 


0.1789 ±0.0080 


-44.0 ± 2.4 


0.59 


891.66 


50.8 


BW 


Lj(782) 


0.0391 ±0.0016 


114.8 ±2.5 


0.56 


782.6 


8.5 


BW 


/o(980) 


0.469 ±0.011 


213.4 ±2.2 


5.81 


975 


44 


BW 


/o(1370) 


2.32 ±0.31 


114.1 ±4.4 


3.39 


1434 


173 


BW 


/2(1270) 


0.915 ±0.041 


-22.0 ± 2.9 


2.95 


1275.4 


185.1 


BW 


Ko*(1430) 


2.454 ± 0.074 


-7.9 ±2.0 


8.37 


1412 


294 


BW 


K^{U30) DCS 


0.350 ± 0.069 


-344. ± 10. 


0.60 


1412 


294 


BW 


K2*(1430) 


1.045 ± 0.045 


-53.1 ±2.6 


2.70 


1425.6 


98.5 


BW 


^1(1430) DCS 


0.074 ± 0.038 


-98 ± 30 


0.01 


1425.6 


98.5 


BW 


K*(1410) 


0.524 ±0.073 


-157 ± 10 


0.39 


1414 


232 


BW 


K*(1680) 


0.99 ±0.31 


-144 ± 18 


0.35 


1717 


322 


BW 


p(1450) 


0.554 ±0.097 


35± 12. 


0.28 


1406 


455 


GS 


0"! 


1.346 ± 0.044 


-177.5 ±2.5 


9.11 


484 ±9 


383 ± 14 


BW 


0-2 


0.292 ±0.025 


-206.8 ±4.3 


0.98 


1014 ±7 


88± 13 


BW 


Non resonant 


3.41 ± 0.48 


-233.9 ±5.0 


6.82 









Table 1: Amplitudes, phases and fit fractions of the different components obtained from the 
likelihood fit of the — > i^s-7r~7r~^ Dalitz distribution in D*^ — > D^tt^ data. Masses and widths 
of all resonances except cJi and o"2 are taken from The abbreviations BW and GS stand for 
relativistic Breit-Wigner and Gounaris-Sakurai ^Hl respectively. The total fit fraction is 1.24. 
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D° mass (GeV/c^) 



Figure 2: The Kstt tt^ invariant mass distribution in the flavor tagged D* sample. 



regions of interference between doubly Cabibbo suppressed and Cabibbo allowed decays and CP 
eigenstate decays exhibit the highest sensitivity to 7. 



5 CP FIT TO B ^ D^*^^K SAMPLES 



A maximum likelihood fit {CP fit) is performed on the B D^*^'^K samples to extract si- 



multaneously the CP violation parameters 7, 6^^\ and r^^' and the signal and background yields. 
The likelihood for each candidate j is obtained by summing the product of the event yield Ni and 
the probability Vi over the signal and the three background hypotheses. The extended likelihood 
function is 



(*) 



c 



exp 



(4) 



The probabilities Vi are evaluated as the product of the PDFs for each of the independent variables 
Xj = {r^Es, AE, J^}. "PP^^ (m^ , ) is the PDF for the Dalitz distribution for the i*^ category. The 
categories in the fit are signal B~ D^*^'^K~ , the continuum background, BB background, and 
B~ — > D*^7:~ and are shown in Fig. ^ The tues and AE distributions for signal events are 
described by a Gaussian. The Fisher PDF is parametrized with a double Gaussian function. The 
signal PDF parameters are determined from the B~ D^*^^7r^ control sample. 



5.1 Background Composition 

The numbers of events for the various background components in the B — > D^*^^K^ samples are 
summarized in Table |21 The dominant background contribution is from the random combination 
of a real or fake D^*)" meson with a charged track in continuum events or other BB decays. 
The combinatorial background in the mES distribution is described by a threshold function [14| 
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Figure 3: (a) The ^ Kstt~^tt~ Dalitz distribution from the D*~^ — ^ D^tt'^ events. Projections 
on (b) m\, (c) m?_, and (d) m'^{Tr^7r~) are shown. The result of the fit is superimposed as a sohd 
line. 



whose free parameter ^ is determined from the B~ D^*^^tt^ data sample. The shape of the 
combinatorial background niES distribution in generic BB decays is taken from simulated events. 
The AE distribution is described by a straight line whose slope is extracted from a fit to the 
B~ D^*^^Tr~ sample. The PDF of the Fisher distribution for continuum events is determined 
from the hies sideband in the same data sample. The Fisher PDF for BB events is assumed to be 
the same as that for the signal. 

An important class of background events arises from continuum where a real D° is produced 
back-to-back with a kaon. Depending on the flavor-charge correlation this background can mimic 
either the b ^ c or the b ^ u signal component. In the likelihood function we take this effect 
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m,2(GeV-/c'*) 



Figure 4: L) — > Kst^^'k Dalitz distribution of simulated {B^ D K^) events. Each event is 
given a weight ^^r^- The black points represent the same events with weight equal to unity. 



into account with two parameters, the fraction fj^o of background events with a real D and the 
parameter R, the fraction of background events with a real associated with an oppositely 
flavored kaon (same charge correlation as the b ^ u signal component). 

The fraction of real from continuum events has been evaluated in events satisfying Tn.ES 
< 5.272 GeV/c^ after removing the requirement on the mass. The fraction R of background 
events with a genuine associated with a negatively charged kaon is obtained from simulated 
events. The values of f£,o and R for continuum qq are summarized in Table Ol The fraction of 
events with a real in generic BB events is found to be few percent. 

A small background originates from B~ — > D^*^^-k~ where the bachelor pion is misidentified 
as a kaon. These events have the same niES distribution as signal but can be distinguished using 
their /S.E information. 

5.2 Likelihood fit on control samples 

We test our CP fit procedure on two high statistics control samples: D*~^ — > D^tt'^ from cc 
continuum events and B~ — > D'^*^^'n~ . The D*~^ — > D^-k'^ sample mimics a B~ D^K~ sample 
with rB= 0. The B~ — > D^7r~ sample is similar to B^ D^K~, but its rs is expected to be 
approximately 0.007 [15,. In the CP fit to D*+ D^ti+ we obtain tb = (-5.2 ± 5.2) x 10"^. 
In the B- D^tt" we obtain rB= (1.8 ± 1.5) x 10"^^ 7 = (18 ± 45)°, 6b = (246 ± 43)° and in 
B- D*'^TT- we find r%= (4.6 ± 2.1) x 10"^, 7 = (90 ± 35)°, S% = (117 ± 35)°. The results 
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Background components 


D^K- 


(D'Jvr^^) K- 




Continuum qq 


125 ±6 


9±2 


38 ±3 


BB 


28 ± 7 


4±2 


14 ±5 


Dtt 


9±8 


0±5 


6±4 



Table 2: Estimates of the numbers of background events in the rriEs > 5.272 GeV/c^ region 
from the fit to data in the fuh ruES region for the B' D^K~, B~ D*^{D'^tt^)K~ , and 
B- D*°{D°j)K- samples. 



qq background parameters 


D^K- 




{D^-f)K- 




0.27 ±0.06 


0.27 ±0.13 


0.13 ±0.02 


R 


0.21 ±0.03 


0.23 ±0.11 


0.16 ±0.06 



Table 3: qq background parameters for the B D^K , B — > D*^{D^tt^)K , and B — > 
D*°{D^-/)K- samples. 

obtained are consistent with the expectations of Monte Carlo experiments. 

5.3 Likelihood fit on D^*^^K- sample 

In the sample of 227 million BB events we obtain the following signal yields 

N{B- ^ D^K-) = 282 ±20, (5) 
N{B- ^ D*^(D^i:^)K-) = 89 ± 11, 
N{B- ^ D*^{D^-l)K-) = 44 ±8, 

in agreement with our expectation from simulation and measured branching ratios. We obtain the 
following CP parameters from the fit for B' D°K~, rB= 0.117 ± 0.053, 6b = (109 ± 28)°, and 
7 = (66 ± 28)° and for B- D*^K-, r%= 0.167 ± 0.065, 6% = (294 ± 28)°, and 7 = (68 ± 29)°. 
These errors are estimated with a Gaussian assumption for the likelihood. However, for this small 
sample, these low rs and r*^ fitted values lead to a non-Gaussian behavior of the likelihood 
function as shown in Fig. El necessitating a different approach, described next, in the computation 
of the confidence intervals for r^, 7 and 5b- Fig. Elshows for tb values generated in the [0.0,0.3] 
range the values obtained in fits to Monte Carlo experiments of the same size as data. While 
the CP fit is linear for large values of r^, it is not sensitive to values below 0.1. This problem 
did not exist for the larger D*'^ — > D^Tr"*" and B~ D^*^^Tr~ samples as we have verified using 
Monte Carlo simulation. 

In Fig. [3 and Fig. |H1 we show the Dalitz distribution and the and m?_ projections for 
events with uies > 5.272 GeV/c^ for B~ D^K~ and B~ -> D*^K~ respectively. S+ and B~ 
candidates distributions are separately shown with the total PDF superimposed. 

5.4 Confidence intervals of the CP parameters 

We evaluate the likelihood function C (r^, 7, 5b) after fixing all parameters except (r^, 7, 5b) 
which are varied in their range of definition, [0,1], [-vr,7r], and [0,27r], respectively. We then estimate 
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Figure 5: InC contour plots in 7 versus in — > D^K (top) and B — > D*^K (bottom). 
Sb is fixed to the value obtained from the fit. Each contour represent a ln>C variation of 0.5. 
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Figure 6: rs values obtained in the likelihood fit versus generated rs values in Monte Carlo 
D^K~ experiments. The error bars represent the RMS of the values returned by the 
fit. The dashed horizontal line indicates the rs value (0.106) found in data. 



the confidence region for the CP parameters using a Bayesian technique. This implies a choice of a 
priori distribution. For this preliminary result we arbitrarily assume a uniform a priori distribution 
for each of the CP parameters ^^,7 and 5b- 

In 'y-TB space we define a two-dimensional confidence region T>[C) corresponding to a given 
confidence level C 

Id(C) '^"'b'^'iC dSBC{rB,y,SB) ^ ^ 

We uniquely define D(C) by requiring that the likelihood value at any point on the boundary of P 
be the same and integrating over all likelihood values larger than the value at the boundary. 

Similarly we define one dimensional confidence intervals 1(C) corresponding to a confidence level 
C. For example, the interval for rs is defined as 

The two dimensional confidence regions I'(C) in 7 versus and 7 versus are shown in 
Fig. 1^1 The red (dark) and yellow (light) regions correspond to the 68% and 95% confidence levels 
respectively. The likelihood distributions for r^, 7, and 6b obtained by integrating /l{rB,^, Sb) 
over the other two variables are shown in Figs. IIUI 1111 and respectively. The 68% and 95% 
confidence intervals are shown in red (dark) and yellow (light) respectively. The intervals for 7 and 
6b are disjoint as a consequence of the 7 ^ 7 it tt and 6b ^ ± vr ambiguities. 

5.5 Constraints on tb, 5b and 7 

From the procedure above we obtain Bayesian 68% confidence intervals. We quote as the central 
values for 7 and 6b the average values weighted by the likelihood distribution. The errors associated 
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Figure 7: Dalitz distribution (first row) from the B' D^K' events with tues > 5.272 GeV/c^. 
Projections on m\ (second row) and m?_ (third row) are shown with the result of the fit superim- 
posed. In the left column candidates are shown, in the right B~ candidates. 
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Figure 8: Dalitz distribution (first row) from the B~ — > D*^K~ events with niES > 5.272 GeV/c^. 
Projections on (second row) and m?_ (third row) are shown with the result of the fit superim- 
posed. In the left column candidates are shown, in the right B~ candidates. 
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Figure 9: Bayesian confidence regions for 7 versus rs (top) and r*^ (bottom) for the D^K~ 
and D*^K~ samples respectively. The red (dark) region corresponds to the 68% confidence 

level region, the yellow (light) region to the 95% C.L. region 




Figure 10: Probability density function for (left) and r*^ (right). The red (dark) region 
corresponds to the Bayesian 68% confidence level region, the yellow (light) region to the 95% C.L. 
region. 
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Figure 11: Probability density function for 
sample . The red (dark) region corresponds 
yellow (light) region to the 95% C.L. region, 
the 7 ^ 7 lb vr and Sb — > 5^ ± vr ambiguities. 




7 in B- D^K- (left) and D*^K- (right) 

to the Bayesian 68% confidence interval region, the 
The intervals for 7 are disjoint as a consequence of 




Figure 12: Probability density function for 5b in B~ D^K~ (left) and B~ D*^K~ (right) 
sample. The red (dark) region corresponds to the Bayesian 68% confidence interval region, the 
yellow (light) region to the 95% C.L. region. The intervals for 6b are disjoint as a consequence of 
the 7 ^ 7 lb vr and 5b — > (5b ± vr ambiguities. 
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with the central values are defined by the boundaries of the 68% confidence interval. We obtain 
rB= 0.087±to:o74, 5b = 114° ± 41°(294° ± 41°), 7 = 70° ± 44°(-110° ± 44°) for ^ D^R- and 
r%= 0.155l|]:[]5?, 5|j = 303° ± 34°(123° ± 34°), 7 = 73° ± 35°(-107° ± 35°) for ^ D*^K-. We 
constrain to be < 0.16 at 90% confidence level. 

As illustrated in Figs. ISland lT^ the data have no sensitivity to 7 and 6b for small values of r^. 
Those values are not excluded by our data. 

We construct a combined likelihood function from the product of the individual likelihoods for 
B~ — > D^K~ and B~ D*^K~ and we repeat the procedure outlined above. From this we obtain 
7 = 70° ± 26° (-110° ± 26°). Fig. ITl shows the 7 likelihood distribution. 




p r , , I , I I , , , I I I , , I , I I , I , , I I I , , , I I , i 
0.05 0.1 0.15 0.2 0.25 0.3 

Figure 13: Statistical uncertainty in Gaussian hypothesis for 7 as a function of the value of 
VB obtained by the likelihood fit in Monte Carlo experiments. 



6 SYSTEMATIC UNCERTAINTIES 

The principal systematic uncertainty on the measurement of 7 comes from the choice of the model 
used to describe —>■ Kstt^tt^ decay. We evaluate this uncertainty by considering alternative 
models. For each model we generate events and fit both the alternative model and the nominal 
model (defined in Section^ to these events. We quantify this uncertainty using the differences in 
the fitted values for r^, 5b and 7. For models where the /9(1450), the -fir*(1680) and/or the doubly 
Cabibbo suppressed i^Q(1430) and -fC|(1430) are removed or a different description of resonances 
is used, the of fit is not significantly different from that of the nominal model. For these 
models the biases on 7 and vb are negligible and the RMS of the distribution of the differences is 
at most 1° and 0.002 for 7 and rs respectively. As an extreme we consider a model without the 
o"! and/or (T2 scalar, or the CLEO Model ^5]. Fits to these models result in a significantly larger 
than that of the nominal model. To illustrate the magnitude of this variation Fig. 1151 shows the 
result of a fit with the CLEO model. For these extreme models the biases for CP parameters are 
still small. The RMS of the differences for 7 and vb are approximately 10° and 0.02 respectively. 
We conservatively assign a.y,asg = 10° and = 0.02 as systematic uncertainties associated with 
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Figure 14: Probability density function for 7 from the combined samples. The red (dark) region 
corresponds to the Bayesian 68% confidence interval region, the yellow (light) region to the 95% 
C.L. region. 



the Dalitz model. 

The summary of the estimates of other systematic uncertainties is given in Table |3J The 
most important effect is due to the uncertainties on the knowledge of the Dalitz distribution of 
background events and of the mES-, and J- PDF parameters for both background and signal. 
Uncertainties in the efficiency variation across the Dalitz distribution are estimated. The statistical 
uncertainty on the amplitude and phases of the nominal Dalitz model also contributes significantly 
to the systematic uncertainties on the CP parameters. 





B- D^K- 


B- D*^K- 


Source 




7 




^B 


7 




Combinatorial background Dalitz shape 


0.008 


6.7° 


3.3° 


0.010 


2.9° 


5.1° 


mESi ^E-, ^ PDF shapes 


0.007 


5.4° 


4.2° 


0.025 


1.8° 


8.2° 


R 


0.018 


3.1° 


3.0° 


0.018 


3.1° 


3.0° 


Efficiency 


0.004 


3.0° 


2.7° 


0.005 


3.0° 


2.8° 


Dalitz amplitude and phase uncertainties 


0.004 


1.6° 


4.7° 


0.014 


6.1° 


8.8° 


Total 


0.022 


9.8° 


8.3° 


0.036 


8.2° 


13.8° 



Table 4: Summary of the contributions to the systematic errors on rs, 7 and 6b- 



7 RESULTS AND SUMMARY 

We report preliminary results of the measurement of and of the angle 7 using the B~ meson 
decays into D^K^ and D*^K^ with a technique based on the Dalitz analysis of the Kstt^tt^ 
three-body decay. From 227 million BB pairs collected by the BABAR detector, we reconstruct 
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Figure 15: CLEO model fit: a) the D Kstt^tt" Dalitz distribution from the D*~ D it" . 
Projection on (b) rn^, (c) and (d) rn^(7r"'"7r~) are shown. The result of the fit is superimposed 
and it is used to estimate the uncertainty due to the Dalitz model. 



282 ±20 B- D^K-, 89±11 5" ^ D*^K-, D*° ^ and 44±8 ^ D*^K-, D*° ^ D°7 

signal events. 

Values of the ratio oi b ^ u and b —>^ c amplitudes for the processes D^K~ and B~ 

D*^K~ at the small end of our measurements allow no determination of 7 at this statistical level. 
Accounting for systematic uncertainties, we constrain these ratios to be tb < 0.19 at 90% confidence 
level and r*g = ± 0.040 lb 0.020. The relative phases between these two amplitudes are 

5b = 114°±41°±8°±10°(294°±41°±8°±10°) and 6% = 303°±34°±14°±10°(123°±34°±14°±10°). 
The first error is statistical, the second error accounts for experimental uncertainties and the third 
error reflects the Dalitz model uncertainty. By combining the information from the two samples 
we obtain 7 = 70° ± 26° ± 10° ± 10° (-110° ± 26° ± 10° ± 10°). For this preliminary result we have 
quoted confidence intervals obtained with a Bayesian technique assuming a uniform prior in rs, 7 
and Sb- 
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